function out = f_angle_class_center(C,F_centers,projection1)
F_centers
if projection1==1
    for i=1:25
%         Index = find(save_data(:,3)==i);
        center_clean(i,:) = F_centers(i,:);
        center_clean(i,:) =center_clean(i,:)./sqrt(center_clean(i,1)^2 + center_clean(i,2)^2); % projecting
    end
    C = C./sqrt(C(1)^2 + C(2)^2);
else
    for i=1:25
%         Index = find(save_data(:,3)==i);
        center_clean(i,:) =  F_centers(i,:);
    end
end

for i=1:25
    dist(i) = sqrt(   (C(1)-center_clean(i,1)).^2 +   (C(2)-center_clean(i,2)).^2 );
end
[Y minIndex] = min(dist);
out.class = minIndex;
